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ABSTRACT 

Star formation from matter including the hot CNO processed ejecta of asymptotic 
giant branch (AGB) winds is regarded as a plausible scenario to account for the 
Q I chemical composition of a stellar second generation (SG) in Globular Clusters (GCs). 

H ' The chemical evolution models, based on this hypothesis, so far included only the 

C/3 ' yields available for the massive AGB stars, while the possible role of super-AGB 

^, I ejecta was either extrapolated or not considered. In this work, we explore in detail the 

role of super-AGB ejecta on the formation of the SG abundance patterns using yields 
recently calculated by Ventura and D'Antona. 
^ , An application of the model to clusters showing extended Na~0 anticorrelations, 

CS| ' like NGC 2808, indicates that a SG formation history similar to that outlined in our 

0^ I previous work is required: formation of an Extreme population with very large helium 

O^ . content from the pure ejecta of super-AGB stars, followed by formation of an Inter- 

\l^ ' mediate population by dilution of massive AGB ejecta with pristine gas. The present 

ff^ I models are able to account for the very 0-poor Na-rich Extreme stars once deep- mixing 

f^ . is assumed in SG giants forming in a gas with helium abundance Y> 0.34, which 

04 ' significantly reduces the atmospheric oxygen content, while preserving the sodium 

abundance. On the other hand, for clusters showing a mild Q-Na anticorrelation, like 

M 4, the use of the new yields broadens the range of SG formation routes leading to 

abundance patterns consistent with observations. Specifically, our study shows that 

p\^ . a model in which SG stars form only from super-AGB ejecta promptly diluted with 

H ■ pristine gas can reproduce the observed patterns. We bricfiy discuss the variety of 

- - -' (small) helium variations occurring in this model and its relevance for the horizontal 

branch morphology. 

In some of these models, the duration of the SG formation episode can be as 
short as ^ 10 Myr; the formation time of the SG is therefore compatible with the 
survival of a cooling fiow in the cluster core, previous to the explosion of the SG 
core collapse supernovae. We also explore models characterized by the formation of 
multiple populations in individual bursts, each lasting no longer than '^10 Myr. 

Key words: globular clusters:general; stars: chemically peculiar; stars:abundances 



1 INTRODUCTION ICarretta etail (12009^ 1 of about two thousand stars in 19 

GCs. Each cluster displays a sodium - oxygen anticorre- 
The presence or multiple stellar populations, in all elob- i j_- ^^i i ^i ^ • r ii j_- i j_- i-rr 

,„„ ^ r 1 lation (though the extension ot the anticorrelation diiter 

ular clusters (GCs) so tar spectroscopically examined, is „ i . . i . \ i , .i . . r i 

^ ' , ,r trom cluster to cluster) due to the existence ol a popula- 

shown with convmcine evidence m the recent analysis oi , . ,. . i • i • i ■ i • i 

tion ol stars which are richer m sodium and poorer m oxy- 
gen than halo stars with the same metallicity. The Na-O 
* E-mail: annibale.dercole@oabo.inaf.it anti-correlation is typical of GCs, whose constituent stars 

© 2011 RAS 



A. D'Ercole et al. 



belong to two or more stellar populations differing in the 
abundances of the elements produced by the hot CNO cy- 
cle and by other proton-capture reactions on light nuclei. 
In fact, these chemical signatures are present also in turn- 
off stars and among the subgiants (e.g. iGratton et al.ll200ll : 
iBrilev et aDl2002l . [2004h . so they can not be due to "in situ" 
mixing in the stars, but must be due to some process of 
self-enrichment occurring at the first stages of the cluster 
life. 

Photometric evidences for the presence of multiple pop- 
ulations are also numerous, and sometimes suggestive of star 
formation occurring in separate successive bursts. The pho- 
tometric signatures of different populations can be ascribed 
in part to helium differences, inferr ed from the morphol- 
ogy of the horizontal branches (HB) |D'Antona et al.l 12002 : 
iD'Antona fc Caloill2004l : iLee et allbOOSl ), or from th e pres- 
ence of multiple main sequences (MS), in m Cen (INorris 
_2004; Piotto ct al. 200^) and NGC 2808 (ID'Antona et al 
20051 : iPiotto et al.ll2007l '). The strong link between the abun- 



dance anomalies (sodium rich and s odium poor groups) in 
the cluster M 4 IIMarino et al.l 20081) an d their location on 
the HE is proven bv lMarino et alj 1 20111 ). who find that the 
blue HE stars of M 4 have high sodium, while the red HE 
stars have normal sodium. This supports the interpretation 
that there is a slight increase in the helium abundance of the 
"ano malous" stars (the sodium rich group) (jD'Antona et al.l 
|2002| ). 

Two other important links between photometric and 
spectroscopic evidence have been recently given. The first 
one comes from the analysis of elemental ab undances of two 
main sequence (MS) stars in NGC 2808 (jBragaglia et al.l 
l20ld ) showing that the blue MS star only has "anoma- 
lous" sodium and aluminum. This supports the interpre- 
tation that the blue MS was formed from helium rich, CNO 
processe d gas. The second li nk has been stressed most re- 
cently bv lPasguini et al.l (|2011i ) who, directly comparing the 
He I 10 830 A lines in the spectra of two NGC 2808 gi- 
ants having different O and Na abundances, find direct ev- 
idence of a significant He line strength difference. From a 
detailed chromospheric modeling, they show that the differ- 
ence in the spectra is consistent and most closely explained 
by an He abundance difference between the two stars of AY 
>0.17, consistent with the expected difference in abundance 
required by stellar models to account for the blue MS. 

Based on these and on a large body of other complex 
evidences, the formation of globular clusters is now schema- 
tized as a two-step process, lasting no longer than ~100 
Myr, during which the nuclearly processed matter from a 
"first generation" (FG) of stars gives birth, in the cluster 
innermost regions, to a "second generation" (SG) of stars 
with the characteristic signature of a distribution of element 
abundances fully CNO-cycled. 

Two main scenarios have been proposed to identify the 
matter constituting the SG stars: these could have formed 
either from the gaseous ejecta of massive asymptotic gi- 
ant branch (AGE) stars ("AGE scenario" : IVentura et al.l 
I2OOII : ID'Antona fc Caloill2004l : iKarakas et al.ll2006l ') or fi-om 
the eje cta of fast rotating massive stars ( "FRMS sce- 
nario" : iPrantzos fc Charbonnell l2006l : iMevnet et all I2OO6I : 
iDecressin et al.l 12007 ) . These studies investigated whether 
the observed abundance patterns found in GCs can be ac- 



counted for by models. Problems are present in both scenar- 
ios (see, e.g.,|Reirzirii 2008, for a critical review). 

In order to go from a "scenario" to a "model", the SG 
formation must be investigated quantitatively, and must rely 
as far as possible on computed models. From this point 
of view, very few are the "models" available in the litera- 
ture, in spite of the large number of works recently pub- 
lished on this subject. In particular, the O-Na anticorrela- 
tion is a signature present in all cluster, but the cluster-to- 
cluster differences are very large: some clusters (the most 
massive ones) show very extended anticorrelations, down to 
values of [O/Fe] , more than 1 dex smaller than the values of 
the FG stars, other clusters show a "mild" anticorrelation, 
in which [O/F e ] is re duced by ~0.2-0.3 dex. In addition, 
ICarretta et al.l (|2009al ) point out that there is a direct cor- 
relation between the minimum [O/Fe] and the maximum 
[Na/Fe] for the examined clusters. Any model, in the end, 
mus t account quan t itativ ely for these features. 

iD'Ercole et al.l l|2008l l followed the hydrodynamical for- 
mation of a cooling flow at the cluster center that occurs 
following the epoch of supernovae type II (SN II) explosions 
LJ, and the associated loss of the remnant gas from which 
the FG stars had formed. The cooling flow is due to the 
low-velocity stellar winds and planetary nebulae of super- 
AGE and massive AGE stars, and meets the physical condi- 
tions for a second epoch of star formation. They show that 
the three populations with different helium content, neces- 
sary to explain th e presence of a triple main sequ ence in 
NGC 2808 (D' Antona et al.ll2005l : IPiotto et al.ll2007l ) can be 
reproduced if the most helium rich population forms from 
the pure ejecta of super-AGE stars, collecting in the cluster 
core devoid of pristine gas after the end of the SN II epoch. 
Afterwards, pristine gas is re-accreted and mixes with the 
massive AGE ejecta, giving origin to the Intermediate pop- 
ulation having a helium content intermediate between that 
of the AGE and of the pristine gas. This hydrodynamical 
model could investigate only the helium content behaviour, 
as the yields of super-AGE stars for elements different from 
helium were not available at that time. 

It is important to understand whether the AGE pollu- 
tion model is c onsistent with the oth er abundance patterns 
of GC stars, so ID'Ercole et all (|20ld ') (hereafter D2010) de- 
veloped a chemical evolution model to test how the O-Na 
and Mg-Al anticorrelations c ould be built up in d i fferent 
clusters. D2010 adopted the I Ventura fc D'Antonal (|2009l ') 
yields for massive AGE ejecta, empirically modified and ex- 
trapolated by educated guesses to account for the super- 
AGE ejecta. D2010 reproduced not only the abundance pat- 
terns of NGC 2808, but also the very different O-Na abun- 
dance pattern of a cluster with a mild O-Na anticorrelation 
like M 4. The conditions required to reproduce the mild O- 
Na anticorrelations are : 1) that the Exreme population from 
pure AGE ejecta could not form in this cluster (and/or in 
clusters with similar patterns); 2) that the dilution with pris- 
tine gas, and the consequent intense SG star formation due 
to the increased amount of available gas, had to be delayed 
by ~30 Myr after the end of the SN II epoch, a constraint 



^ Actually, we refer here to the explosion of core-collapse super- 
novae, thus including also type lb and Ic; however, for brevity, in 
the rest of the paper we refer to all these as SN II. 
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that should be modeled and further explored with appropri- 
ate dynamical simulations. 

In this paper we further extend the models presented in 
our previous work to include the new yields for super-AGB 
model s from 6.5 to 8 Mq computed bv lVentura fc D'Antonal 
l|2011l ). In Sect. [2] we outline the new complete set of yields 
now available and introduce a hypothesis of in situ "deep- 
mixing", in the high helium red giants, that helps to deal 
with the extreme oxygen anomalies. In Sect. [3] we summa- 
rize the D2010 framework and its relevant parameters. We 
present our models for the chemical evolution of M 4 in 
Sect. H and of NGC 2808 in Sect. \5\ Interestingly, models 
based on the new extended set of yields lend further sup- 
port to our previous models while solving some of their dif- 
ficulties. In some of M 4 models the second star formation 
epoch may occur in just ~10 Myr, so these models are not 
in contrast with the formation of SG massive stars that, ex- 
ploding as SN II, prevent any further star formation. We 
discuss in Sect. l5.3l models including an SN II phase also for 
the SGs of NGC 2808. We summarize the results and discuss 
the implications for the initial mass of the FG in Sect (6] 



2 SUPER-AGB AND AGB YIELDS 

In Table [T] we list the average abundances in the 
ejecta of AGB and super-AGB stars adopted in 
this work. These yields have been calculated by 
I Ventura fc D'Antonal (|2009l ') (see their Table 2) and 
further extended with the calculations of t h e yiel ds for 
super-AGB stars by IVentura fc D'Antonal (|201ll ') and 
I Ventura. Carini fc D'Antonal ( 20111 ) . The latter works show 
the continuity between the two set of yields. The models 
have been computed by using the same input physics 
and nuclear reaction rates, so the results can b e matched 
together. We refer to lVentura fc D'Antonal (|201ll ) also for a 
comparison with other super-AGB results in the literature 
(notably bv iSiesjIJOlOJ '). Fig. [l] shows the O-Na yields for 
these models, together with the yields adopted in D2010. 

Notice that the maximum mass that does not ignite 
as SN II in the super-AGB models here adopted is 8 A'Iq . 
Thus the second generation formation in these new chemical 
evolution models can not start before ;$ 40Myr from the 
FG birth, as this is the evolutionary time computed for the 
8 Mq star. In D2010 we assumed that super-AGB would 
include the 9 Mq star, evolving at ~32 Myr. The minimum 
mass leading to SN II explosion (and thus the duration of 
the SN II epoch) depends mostly on the assumptions made 
on the core overshooting during the core hydrogen burning. 
Our models employ a moderately large overshooting (e.g. 
ISchaller et aLlll992l : IVentura et al.lll998l '). If this is reduced, 
the minimum mass of possible SN II progenitors increases, 
and the lifetime of the SN II epoch decreases. This input 
defines not only the beginning of the phase during which the 
low velocity winds accumulate for the SG formation (after 
the FG SN II stage), but also the duration of a possible, 
subsequent, SG SN II phase (from ~30 to ^40 Myr, see 
Sect. [Oil. 
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Figure 1. The yields of metallicity Z=10~ from Ventura & 
D'Antona (2009) from AGB models are plotted as full black 
squares. Masses from top to bottom are 3, 3.5, 4, 4.5, 5, 5.5, 6 and 
6.3 Mq. Yields from Ventura & D'Antona (2011) for super-AGB 
models of 6.5, 7, 7.5 and 8 Mq are plotted as open black squares. 
Red triangles represent the yields adopted in D'Ercole et al. 2010: 
masses are 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.3, 6.5 and 9Mq. 



2.1 Sodium 

By increasing the initial mass from 3 Mq to 8 Mq , the Na 
and O values provided by the models ejecta first decrease, 
reach a minimum value in the yields (at masses ^6.5 Mq), 
and then they both revert back toward higher values. No- 
tice that the location in the Na-O plane of the yields of 
super- AGBs for the masses 6.5^M/Mq^8 is not different 
from that of stars with lower masses, it is simply reversed. 
This b ehaviour is due to the huge mass loss rates of super- 
AGBs (|Ventura fc D'Antonall201ll ) and differs from the val- 
ues adopted in D2010 in which the last point in the yield 
table was taken to be at very low oxygen ([0/Fe]=-l) and 
moderate sodium. 

The comparison between yields and the data for single- 
metallicity GCs clearly shows that to reproduce the observed 
Na-O anticorrelati on it is necessary to dilute the AGB ejecta 
with pristine gas iJBekki et al.l 200?! ; JD'Antona fc Ventura! 



witn prii . 

l2007l : (D'Ercole et al.ll2008l . l201ol ) 

Although the yields of massive stars do not show 
the strong direct correlation of Na and O, but a mod- 
erate anticorrelation, diluti on is also needed in models 
based on the FRMS electa JPrantzos fc Charbonnell 120061 : 
JDecressin et al.l l2007l : JLind et al.l |201l|), and seems to be 
present in all clusters. Detailed hydrodynamical models able 
to explain the dynamics and the presence of pristine gas dur- 
ing the SG formation ar e still lacking. The iss ue of dilution 
is discussed in detail bv JD'Ercole et al.l (|201ll ). 

The Na abundances adopted in D2010 are -1-0.3 dex 
larger than the values reported in IVentura fc D'Antonal 
(2009) and adopted in the present work. This assumption 
was necessary to succesfuUy reproduce the observed anti- 
correlations, and was theoretically justified by assuming a 
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Table 1. Averaged abundances in the ejecta of massive AGB and super— AGB stars. 



M/Mq t/IO^" Mc/Mq^ Y 



[O/Fe] [Na/Fe] [Mg/Fe] [Al/Fe] loge(Li)'= 



3.5 

4.0 

4.5 

5.0 

5.5 

6.0 

6.3 

6.52 

7.0 

7.5 

8.0 



332 

229 

169.5 

130.3 

103.8 

85.1 

71.2 

65.2 

61.5 

53.7 

46.8 

38.8 



0.76 
0.80 
0.83 
0.86 
0.89 
0.94 
1.00 
1.03 
1.08 
1.20 
1.27 
1.36 



0.248 
0.265 
0.281 
0.310 
0.324 
0.334 
0.343 
0.348 
0.352 
0.358 
0.359 
0.344 



0.92 

0.77 

0.44 

0.19 

-0.06 

-0.35 

-0.40 

-0.37 

-0.24 

-0.15 

0.01 

0.20 



1.16 
1.30 
1.18 
0.97 
0.60 
0.37 
0.31 
0.30 
0.32 
0.39 
0.67 
1.00 



0.57 
0.55 
0.48 
0.43 
0.35 
0.28 
0.29 
0.30 
0.34 
0.36 
0.415 
0.440 



0.65 
0.66 
0.55 
0.85 
1.02 
1.10 
1.04 
0.99 
0.91 
0.86 
0.65 
0.50 



2.77 
2.43 
2.20 
2.00 
1.98 
1.93 
2.02 
2.06 
2.36 
2.12 
2.75 
4.39 



° Total evolutionary time until the AGB phase 
''Core mass at the beginning of the AGB phase 
'=loge(Li)=log(ArLi/iVH) + 12. 

1 Yields for 3s£M/M0s:6.3 from I Ventura 

2 Yields for 6.5^M/Mq^8.0 from I Ventura fc D'Antona (20l3) and I Ventura. Carini fc D'Antonal 1I2OIII ') 
^ The values listed in the present Table are obtained for the following set of initial abundances in mass fraction: H=0.75, He=0.24, 
Li=10-", 0=0.849 X 10"*, N=0.249 x 10""', 0=0.58 x 10-^, Na=0.103 x lO-^, Mg=0.497 x lO"", Al=0.178 x lO'^, Fe=0.377 x 10""'. 



D'Antona 1*200' 



]ll) 



50% increase in the ^''Ne(p,7) reaction rate and a larger 
^'^Ne abundance in the initial gas mixture. In the present 
work, use of super-AGB yields makes this adjustment nec- 
essary only for the specific case of M 4 in which a '^^Ne 
overabundance is actually present (see Sect.|3}- 



2.2 Oxygen 

The computed super-AGB models show that the sodium 
and oxygen yields increase with the mass, and that [O/Fe] 
never reaches values smaller than —0.4 (cf. Table [T] and 
Fig. [1}. This behaviour raises a question on the origin of 
the very low oxygen abundances ([0/Fe]~ —1) found in 
the giants of some clusters (these are the stars classified 
as Extreme second generation by Carretta et al. 2009a). We 
propose to solve the issue with a single additional assump- 
tion. We attribute the very small oxygen abundances to an 
"in situ" deep-mixing, acting only in the progenitors hav- 
ing a very high helium content. In these stars, the small 
molecular weight discontinuity during the red giant branch 
evolution may not be able to preclude deep-mixing and low- 
ering of oxygen in the envelope. Stars having the largest 
helium abundances would then signficantly reduce their at- 
mospheric o xygen content, while preserv ing a similar sodium 
abundance (jD'Antona fc Ventural 120071 11. The presence of 



^ Historically, deep-mixing wras modeled in giants starting their 
life from a standard abundance mixture, and produces a sodium 
increase concomitant to the oxygen depletion, dealing with re- 
gions in which neon has been processed to sodium by proton cap- 
tures. These models were in reasonable agreement with the whole 
sodium versus oxygen anticorrelation (see, e.g. Weiss, Denis- 
senkov & Oharbonnel 2000, their fig. 9), but the deep- mixing 
interpretation had to be abandoned when the chemical anoma- 
lies were shown to be present in turnoff stars. If deep-mixing is 
instead allowed by the high helium content of the Extreme SO 
component, we must take into account that these stars are formed 
from AGB processed matter, in which the initial neon has already 
been converted to sodium. In this case, no additional sodium is 



[O/Fe] abundances smaller than those predicted by AGB 
models only among the giants, a nd not among the tu rnoff 
or subgiant stars so far examined (jCarretta et al.ll2006r ). and 
in stars belonging to clusters also showing the signature of 
a very helium rich MS, such as NGC 2808 and lo Cen, are 
a possible argument in favour of this interpretation. Obvi- 
ously, the determination of a very low oxygen content — 
similar to the one in the most anomalous red giants — in the 
atmosphere of blue MS stars in these clusters would falsify 
this suggestion (jBragaglia et al.ll20ld . were not able to mea- 
sure oxygen in their spectrum of the blue MS in NGC 2808.). 

We apply the following deep-mixing scheme to our mod- 
els: if a population of SG stars is born from the pure ejecta 
of super-AGB stars and massive AGBs, in which the helium 
content is Y>0.34, we assume that the surface abundances 
in the atmospheres of the SG giants is due to deep-mixing; 
as explained above, this has the effect of reducing the oxygen 
abundance while leaving the Na abundance unaltered. 

The following steps are followed: 1) we examine the he- 
lium content in the star-forming gas at each time step; this 
helium content is the average content of the gas reservoir, 
either due to the super-AGB or AGB ejecta, in absence of 
pristine gas, or averaged with the re-accreted pristine gas; 
2) if the helium content in the star forming gas is Y>0.34, 
we adopt a "deep-mixing" scheme to provide the oxygen 
content at the surface of giants. In the absence of non para- 
metric models, the relation chosen to fix the oxygen content 
after the deep-mixing has some degree of arbitrariness. We 
decided to link the value to the helium abundance of the 
gas forming the stars. [O/Fe] is linearly interpolated be- 
tween the values -0.4 and -0.9, assumed respectively for 
helium content Y=0.34 and Y=0.36. Consequently, when 
the ejecta of the mass 7.5 Mq form stars, they will reach 
the minimum oxygen values when evolving as giants. We 
assume: [O/Fe] =-0.4-(Y-0.343)x 37.5. We point out that 



produced inside the star during the h ydrogen burning stages, as 
shown bv lD'Antona fc VenturJ | |2007| ). 
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this expression is adopted as an example only, and other - 
arbitrary- formulations could have been chosen. 3) if the av- 
erage helium abundance is smaller than the values for which 
deep-mixing occurs, the [O/Fe] of the SG stars is equal to 
the average between the standard [O/Fe]~0.4 of the pristine 
matter, and the mildly depleted [O/Fe] of Table [T] 

This unique assumption on the simulations will change 
dramatically our understanding of the O-Na anticorrelation 
in clusters that do not show extreme anomalies. 



2.3 Magnesium and Aluminum 

As for magnesium and aluminum , in this paper 
we adopt the yields obtained in our standard models 
l|Ventura fc D Antonal bOOgj ) in which the NACRE upper 
limits for the al uminum production by proton capture on 
^^Mg and ^''Mg (' Angulo et "allll999l ') are used. The super- 
AGB computations have shown that the maximum Mg 
depletion is achieved in models of 5-6 Mq , and not in 
the super-AGB models, contrary to the extrapolation sug- 
gested in D2010. The extent of possible variation of Mg 
and Al yields by varying the input physics is discussed in 
IVentura. Carini fc D'Antonal (|201ll ). 



2.4 Helium 



As discussed in the literature (|Ventura et al.ll2002l '). the he- 
lium abundance in the ejecta of massive AGB stars is mainly 
due to the second-dredge up, as the thermal pulse phase in 
the mass range of interest for the formation of the second 
generation is so short that the episodes of third dredge up 
do not increase significantly the value achieved at the sec- 
ond dredge up. The same happens in super-AGB stars, 
whose evolution is faster than that of massive AGBs and 
whose thermal pulses are less energetic and less effective 
for what concerns the third dredge up. Therefore, the he- 
lium yield of super-AGB stars, contrary to what happens 
for the other chemical elements of interest here, does not 
depend on the assumptions made for the mass loss rate. 
The average helium abundance increases steadily with the 
mass, reaching a maximum value of mass fraction Y=0.359 
for the 7.5M0model. In the 8 Mq evolution, the total he- 
lium in the ejecta is down to Y=0.344. As discussed in 
IVentura fc D Antonal l|201ll '). this is due to the different be- 
haviour of this model, that ignites carbon in the core before 
the second dredge up. Although the results of the 8 Mq evo- 
lution must be taken with caution for what concerns the 
other elemental abundances, that dramatically depend on 
the huge mass loss suffere d by this mode l (following the 
choice made in [Ventura fc DAntonall201ll . for the formu- 
lation of the mass loss rate), it is difficult to think that a 
decrease of the second dredge up effects is physically moti- 
vated. Therefore, these models can not explain the Extreme 
multiple generations of GCs, if we have to take at face value 
the helium contents Y>0.38 that compete to the blue main 
sequences (|Norrisll20o3 : |Piotto et al.ll2007l 'l or even the value 
Y^O.42 attributed to the blue hook stars in NGC 2419 by 
Idi Criscienzo et al.l l|2011f ). We follow this line of reasoning: 
we must not take at face value neither the fit of observa- 
tio ns with models (see, e.g., the alternative interpretation 
by IPortinari et al.l |201C| . on the Y values to be attributed 



to the helium rich main sequence) nor the mediated propos- 
als of precise very large values of Y from the analysis of the 
horizo ntal branch stars, as often pointed out by the same au- 
thors (jDAntona fc Caloil 120081 ) . but only the existence, in 
these clusters, of well defined helium rich sequences. In the 
absence of direct measurements of the helium abundance, we 
will assume that the average value of Y of our super-AGB 
models is more than adequate to descr ibe the Extr eme pop- 
ulations of these clusters. Notice that ISiesj (|2010|), in spite 
of using different physical inputs, reaches a maximum value 
(Y=0.375) in his super-AGB mode ls of the same metallicity. 
Other models in the literature fe.g. [ Boothrovd fc SackmannI 
1 19991 : lDohertvet"aIll2010l : iLagarde et al.ll2011^ reach values 
of helium abundance in the same range (Y~0.33 - 0.38), 
in spite of differences in codes, numerical inputs and phys- 
ical approaches in the mixing formulation (see for discus- 
sion Ventura 2010 and D Antona, Charbonnel et al. 2012, in 
preparation). 



3 THE CHEMICAL EVOLUTION MODEL 

The chemical evolution is computed following equations 2-5 
of D2010. According to the framework presented in D2010, 
FG stars are already in place and have the same chemical 
abundances of the pristine gas from which they form; the 
SG stars form from super-AGB and AGB ejecta gathering 
in the cluster centre through a cooling fiow and partially 
diluted with pristine gas accreting on the cluster core with 
(possibly) a time delay with respect to the beginning of the 
SG formation epoch. 

In order to preserve the cooling fiow, we still assume 
that only SG stars th at do not explode as SN II can form 
JD'Ercoie et al.ll2008f) FI Another reason to exclude the SG 
SN II is that normally GCs show a very small spread of iron 
(ICarretta et al.ll2009d ). indicating that the SG stars hardly 
have been polluted by SN II ejecta (neither of the FG nor of 
the SG). The parameters characterizing the model are the 
same ones listed in D2010: 

(i) tacc'- time at which the maximum accretion of pristine 
matter occurs; 

(ii) Tpri'- timescale of the pristine gas accretion process; 

(iii) tend: time at which the SG star formation ends; 

(iv) p,,fg: densiy of FG stars; 

(v) po.pr: densiy regulating the amount of available pris- 
tine gas, normalized to p*,fg; 

(vi) u: star formation efficiency 

(vii) X = p.,SG/p*tot: ratio between the SG stars and the 
total nowadays alive stars (p., tot = P*,fg + P*,sg) 

An exploration of this whole set of parameters has been 
presented in D2010. In the following we will vary only the 
most critical parameters (iacc, Tpri, icnd, i^ and po.pr). No- 
tice that fend is not necessarily linked to face and Tpri, as 
the phenomenon causing the end of the SG star formation 
(e.g. the phase of type I supernova explosions, cf. D2010) 
is not related. The ratio x, in most GCs, is larger than 0.5 



^ The following will show that this limitation in the IMF can 
be relaxed for clusters showing a mild O-Na anticorrelation, see 



Sect. [531 



© 2011 RAS, MNRAS OOO.fTHTil 



6 A. D'Ercole et al. 



Table 2. Dynamical input parameters for the simulations 



Cluster 


Model 


SG formation 


tacct") 

(Myr) 


(Myr) 


^end 

(Myr) 


itt) StC'^ 
(Myr) 


P.,FG 
M0pc~^ 


PO.pr 


V 


X 


M4 





continuous 


65 


2 


105 




9.4 


0.091 


0.1 


0.7 


M4 


1 


continuous 


43 


1 


100 




60 


0.085 


0.1 


0.5 


M4 


2 


continuous 


43 


2 


60 




940 


0.05 


0.1 


0.5 


M4 


3 


accum. 


43 


2 


50 


(45) 


940 


0.04 


0.5 


0.5 


M4 


4 


accum. 


43 


2 


55 


(45) 


940 


0.04 


0.5 


0.5 


M4 


5 


accum. 


48 


10 


58 


(48) 


940 


0.08 


0.5 


0.5 


NGC 2808 





continuous 


58 


6 


68 




240 


0.01 




0.5 


NGC 2808 


1 


continuous 


65 


8 


90 




240 


0.0095 




0.4 


NGC 2808 


2 


accum. 


65 


8 


90 


(50) 


240 


0.0095 




0.4 


NGC 2808 


3 


bursts 


87 


10 


100 


50-80 


240 


0.0045 




0.4 


NGC 2808 


4 


bursts 


87 


10 


90 


50-80 


240 


0.0045 




0.4 


NGC 2808 


5 


bursts 


87 


10 


90 


50-80 


240 


0.02 




0.4 



'") times in Myr. 

'*'' if: time until which the AGB ejecta are accumulated before star formation begins 

' ' in models N2808-3, 4 and 5 we list (5t, the time interval during which there is no star formation. 



iJD'Antona fc CaloilliooJ : ICarretta et al.ll2009al ). This cur- 
rent value depends not only on the duration and efficiency 
of the star formation of the SG, but also on the subsequent 
loss of FG stars (see, e.g. iD'Ercole et"ai]|2008l ). Although 
the loss of FG stars does not enter directly into the models 
(it is in fact absent in equations 2-5 of D2010), its value can 
be inferred a posteriori from the fit of specific data for the 
clusters. In the examples of Sect. |4] we adopt x = 0.4 or 0.5 
for NGC 2808 and x = 0.5 or 0.7 for M 4. 

The initial values for the pristine gas abundances can 
be inferred from the observed data. For example, the initial 
Mg and Al abundances in the cluster M 4 must be consis- 
tent with the observed values [Mg/Fe]=0.5 and [Al/Fe]=0.5, 
while we adopt more standard values of [Mg/Fe]=0.3 and 
[Al/Fe]=0 for the initial abundances in NGC 2808. M4 FG 
stars also exhibit a small increase (~ -|-0.05dex) in the ini- 
tial oxygen content. This calls for a variation of the Mg 
and Al abundances of the AGB ejecta given in Table [1] as 
these latter are computed for the initial values [Mg/Fe]=0.4 
and [Al/Fe]=0. Magnesium and Aluminum were modified on 
the basis of the details of the various channels of the Mg-Al 
nucleo synthesis, discussed in I Ventura. Carini fc D'Antonal 
l|201 jl . The change in the oxygen yields is simpler, because 
this element is destroyed with continuity during the whole 
AGB phase, which allows a straight scaling relation. Fur- 
ther, the overabundance of a-elements in M 4 should also 
imply an overabundance of neon. Adopting the relation dis- 
cussed in D2010 for the sodium production by AGB stars 
(S [Na/Fe]~ Jp^Ne/Fe]), we increased the sodium yield val- 
ues of Table [l] only for the specific case of M 4. Table [3] 
summarizes all the variations adopted in our models. 



4 RESULTS: M 4 

4.1 Continuous star formation 

In D2010 we p resented a model to fit the M 4 0-Na and 
Mg-Al data bv lMarino et all (|2008l '). In Fig.[2]we reproduce 
a model with the same dynamics (M 4-0 in Table[2| but with 
the abundances calculated using the new Table[T](and Table 
[3|. The bottom right panel shows the role of ejecta (dash- 



dotted line) and of the pristine gas accretion (long-dashed 
line) and their timing in the detailed formation of the SG 
(dashed line). The bulk of the SG formation occurs when the 
pristine matter is re-accreted, at ~65 Myr, and lasts for ~40 
Myr more. The dashed lines in the left panels show the path 
of the O-Na and Mg-Al abundances along the simulation. 
We see that, when dilution with pristine matter sets in, the 
values tend towards the FG abundances. As we discussed in 
D2010, this choice of input parameters provides agreement 
with the data if the sodium abundance is assumed to be 
larger by -1-0.1 dex, with respect to the yields adopted in that 
work, that is, larger by +0.3 dex with respect to Table [T] 
As discussed in Sect. 12.11 this would mean to assume that 
the neon abundance in the gas forming the FG stars of M 4 
is -1-0.3 dex larger than in other cluster-forming regions, 
in particular in NGC 2808, for which we used directly the 
values in Table (2] The top right panel of Fig. [2] shows the 
helium distribution. The total spread is (5Y~0.04, and could 
be revealed by a careful analysis of the MS. 

The new super- AGB yields broaden the range of possi- 
ble star formation routes that can satisfactorily fit the M 4 
data, showing that different re-accretion histories may lead 
to similar mild O-Na anticorrelations found in numerous 
GCs. In the new models star formation can begin much ear- 
lier: a first example is shown in Fig. [3] (case M 4-1). The 
star formation efficiency in this model has still the low value 
i/=0.1 assumed in the previous case, and no star forms di- 
rectly from super-AGB ejecta (as we see from the distribu- 
tion of the helium abundance). Here (and in all the following 
models of M 4) the adopted excess of sodium abundance is 
reduced to (5[Na/Fe] = 0.2, and we assume 1=0.5 instead of 
0.7 considered in the previous model for the sake of compar- 
ison with the model of D2010. A fraction of 50% of FG stars 
is in fact consistent wi th the fraction of red HB stars given 
in lMarino et al.l (|2008l ). Recent analysis of UV based colors 
of the red giants actually show that the RG branch is split- 
ted into a Na-rich (redder) and Na-normal (bluer) part, in 
balanced numbers (Milone 2011, private communication). 

Model M 4-1 requires the presence of pristine gas since 
the early phases of SG formation to efficiently dilute the 
super-AGB matter. Thus the strong dilution of the high- 
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Table 3. Chemical input parameters for the simulations 



Cluster 


S[0/Fe] 

(a) 


5[Na/Fc] 

{a) 


<5[Mg/Fc] 

(a) 


<5[A1/Fc] 

(a) 


Oi, 

(fc) 


Mgin 


Alin 

(b) 


M 4 
NGC 2808 


0.05 
0.0 


+0.21 
0.0 


+0.03 
-0.05 


+0.05 
-0.10 


580 
580 


60 

40 


6 
2.5 



(1) for all the M 4 models but M 4-0, for which (5[Na/Fe]=0.3. 

(") difference in the AGB yields with respect to Table 1. 

('') mass fraction of the element in the pristine gas, in units of 10~^. 



helium super-AGB ejecta leads to helium abundances not 
much larger than the pristine helium abundance in the gas 
forming the SG. As discussed in Sect l2.2l the oxygen con- 
tent is not reduced by deep-mixing: the SG ends up forming 
a group of stars with relatively large sodium (also thanks 
to the inclusion of (5[Na/Fe]=+0.2 in the yield table, as dis- 
cussed above), but with oxygen depleted only by ~0.1 dex. 
The small helium spread shown by this simulation is prob- 
ably undetectable on the main sequence, but can lead to a 
distribution of the HB stars such that the Na -rich stars are 
bluer than the Na-normal stars, as found bv lMarino et al.l 
l|201ll ). 

We also run the model M 4-2 (Fig. |4]) to illustrate how 
acceptable solutions can be obtained even for significantly 
shorter evolutionary times. This shows that this kind of 
models are flexible enough to reproduce the observations 
in a variety of conditions; this also means, however, that 
these models can provide only limited constraints on the 
gas hydrodynamics and the related star formation history. 




0.2 0.4 0.6 

[Mg/Fe] 








0.2 0.4 0.6 

[0/Fe] 



40 60 80 100 
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4.2 Star formation after gas accumulation 

In Fig. [S] we model a different case (M 4-3). Here we adopt 
1^=0.5, and let the super-AGB ejecta and pristine gas accu- 
mulate up to an age of 45 Myr, when a burst of formation 
occurs, lasting until a total age fend=50 Myr. The process of 
SG formation must stop very soon with these parameters. 

In Fig|S] (case M 4-4), we assume iond=55 Myr, and 
a small tail of stars with larger helium content (and lower 
[0/Fe]) is formed from pure AGB ejecta in the time interval 
50-55 Myr; such a tail is not present in the data, but it 
is possible to increase Tpri to avoid this problem. Our aim 
here is to show how much a small variation of the input 
parameters may lead to different results, when SG stars form 
from super-AGB ejecta, i.e. when, contrary to model M 4-0, 
shorter evolutive times are involved. 

To further illustrate this point, we show another sim- 
ulation in Fig. [7| (Case M 4-5), with a longer timescale for 
the SG formation. Interesting to notice, the helium content 
of the SG is again only slightly increased, but here it is well 
separated from the FG value. Possibly, a careful color analy- 
sis of th e main sequence widt h following the lines described, 
e.g., bv iMilone et al.l (|2011ah could reveal this bimodality 
in helium abundance, and discriminate between this model 
and the previous one. Models similar to this one may be 
able to reproduce the abundances in other clusters showing 
a bimodal HB. 



Figure 2. Simulation for the SG of M 4, following the criteria 
of D2010 (case M 4-0). The top-left panel and the bottom-left 
panel show the [Mg/Fe]-[A1/Fe] and the [0/Fe]-[Na/Fe] diagrams, 
respectively. Data (open circles) are from Marino et al. (2008). 
Squares and triangles represent a sampling of the SG and FG 
stars by our model, respectively. The dashed lines represent the 
gas trajectory within the diagrams; the sampled stars in principle 
would be located on this line, but we introduced a random scatter 
in the range 0.1 dex in their coordinates in order to mimic the 
observational errors. The top-right panel illustrates the stellar 
He distribution. In the bottom-right panel the evolution of the 
following quantities is reported: total amount of gas (solid line), 
stellar ejecta (dot-dashed line), pristine gas (long-dashed line), 
SG stars (dashed line). 



4.3 Which model is more adequate? 

Thus, there are two main different modalities in which the 
AGB pollution may lead to a viable model for M 4: either the 
re-accretion of pristine gas is delayed by more than 20 Myr 
(model M 4-0 and D2010), or it occurs early after the SN II 
epoch (models M 4-2 and following). The choice between 
these possibilities can rely upon the distribution of other el- 
emental abundances in the FG and SG. In particular, the 
model M 4-0 shows a small increase in the aluminum abun- 
dance of the SG (Fig.m top left panel) . iMarino et al] (l2008h 
data, and also the recent data bv lViUanova &: Geisleij l|201l[ ). 
do not show a clear aluminum trend. In addition, the D2010 
and M 4-0 models show a finite spread in the helium abun- 
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Figure 3. Model M 4-1. In this model tlie pristine gas accretion 
is greately anticipated with respect to the model M 4-0 of Fig. [2] 
(symbols as in Fig.[2]l. 
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Figure 5. Model M 4-3. The pristine gas and super-AGB ejecta 
accumulate and a short-lasting burst of star formation occurs at 
45 Myr (symbols as in Fig. [2J. 
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Figure 4. Model M 4-2. The strong, early dilution or super-AGB 
ejecta with pristine gas allows the formation of the SO at early 
times aftr the end of the SN II epoch (symbols as in Fig.[2j. The 
time evolution is much shorter compared to models M 4-0 and M 
4-1. 



dance of the two populations, and model M 4-5 a bimodality 
in the helium content. 

The helium abundance may offer another possible way 
to choose among the models, by examining the morphology 
of the HB, and the abundances of sodium a nd oxygen at dif- 
ferent locations along the HB. Originally, iD'Antona et al.l 
l|2002h proposed to explain the color spread in the HB of 
most GCs not as an effect of a mass loss spread, as gener- 
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[0/Fe] 



Figure 6. Model M 4-4. The same input parameters of model 
M 4-3, but the evolution lasts for 5 Myr longer, when the pristine 
gas is already exhausted, and a tail of helium-rich stars is formed 
from pure AGB ejecta (symbols as in Fig. |2). 



ally thought, but as an effect of helium increase from the 
cooler to the hotter side. At a given age, the mass of stars 
evolving off the main sequence is smaller for larger helium 
abundances. This implies that, for a similar mass loss, stars 
with a larger helium abundance end up with a smaller mass 
and thus at a bluer location on the HB. Depending on the 
metallicity, a very small helium increase (from Y=0.24 to 
Y~0.27) is necessary to produce the gap — a real lack of stars 
in the RR Lyr regi on — between the red clu mp and the blue 
HB in NGC 2808 JD'Antona fc Caloill20oA . In clusters like 
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Figure 7. Model M 4-5. The burst of star formation occurs after 
accumulation until an age of 48 Myr: FG and SG stars are born 
with a finite, although very small, difference in helium content 
(symbols as in Fig.[2)l. 
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Figure 8. Sodium versus helium abundance in the labelled simu- 
lations. The line represents the Y-Na temporal path of the model 
yields during the formation of the SG. As for other figures, the 
number of points along different portions of the lines is propor- 
tional to the star formation rate (symbols as in Fig. [2]l. 



M 4, having a more continuous distribution between red and 
blue, with the presence of several RR Lyr stars, the helium 
variations may have been less abrupt. Marino et al. (2011) 
found that the blue HB stars of M 4 have hi gh Na, while the 
red H B stars have normal Na. Furthermore, IVillanova et al.l 
1)20121 ) have recently observed six blue HB stars in M 4, 
finding that they have high Na, low O and helium content 
Y~0.29. A similar location of Na-rich, 0-poor stars at the 
blue side of the HB has been recently discovered by Grat- 
ton et al. (2011) in NGC 2808, confirming the scenario. We 
remark that this same analysis shows a mild anticorrelation 
of the color with the Na/0 ratio within the red clump stars, 
so that, contrary to the simple interpretation, even the red 
side of the HB may be not fully populated by FG stars, but 
also by a few SG stars having a helium content barely larger 
than the FG value. Also in NGC 1851 Gratton et al. (2012) 
have found a few Na-rich stars at the hotter end of the red 
clump. These cases show that extreme care is necessary be- 
fore drawing firm conclusions on the number ratio of FG to 
SG stars from simplified analyses. . 

Figure |S] illustrates the different ways in which helium 
is correlated to sodium in our simulations. While in model 
M4-0 the difference in helium between SG and FG stars is 
5Y> 0.02, other models have a population of SG stars at 
Y~ 0.25 (model M4-2) or 0.255 (model M4-4). Model M 4- 
5 shows that it is possible to obtain a mono-helium SG (in 
our case Y~0.26) with a large sodium spread. From these 
examples, it appears that the models employing the super- 
AGB may be able to explain the presence of SG stars on 
the blue side of the red clump of NGC 185 ]|j (see Sect. 5 for 



We must keep in mind that the case of NGC1851, however, 
is more complex than those we are exploring here. In this clus- 
ter, metallicity and s-process element differences are also present 
and a significant age difference between the two populations is 



the case of NGC 2808). Quantitative simulations of the HB 
in this and other clusters might help to determine which 
of the models presented is the best. While these simula- 
tions are outside the scope of this paper, we point out that 
the models presented can potentially explain a variety of 
different HB morphologies and can be tested by spectro - 
scopic observations s imilar to those bv lMarino et al.l l|2011f ): 
IVillanova et^ ((2012i fl. 

We summarize these results by concluding that: 

(i) The mild 0-Na anticorrelation and the small - if any 
- variation in helium shown by a class of GCs can be due to 
pollution of super-AGB stars into the re-accreted pristine 
gas forming the second generation stars. The needed ingre- 
dient is that the pristine matter falls back to the cluster core 
very soon after it has been totally removed from the cluster, 
during the FG SN II epoch. 

(ii) Another important outcome of the simulations includ- 
ing the super-AGB ejecta is that the SG formation phase 
can have a short duration (~ 10-20 Myr) and does not need 
to extend up to ~ 100 Myr (as found by D2010). The mod- 
els characterized by a shorter SG formation epoch, on the 



necessary to explain the complex features of its CM diagram (see 
Gratton et al. 2012 for further discussion). 

^ We comment, in passing, on the case of NGC 6397, a clus- 
ter having a short HB and a thin MS, both indicative of s carce 
star-to-star differences in helium l ldi Criscienzo et al.ll201Cll ) but 
a significant spread in sodium. This latter was assumed by 
jCarretta et al.ll2010l ) to correspond to the mixture of FG plus 
SG stars, with ~30% of stars belonging to the FG. We can not 
model NGC 6397 with the present yields, as it has a much smaller 
metallicity, but we remark that it can be reanalysed in terms of 
a model similar to the SG of model M 4-5 in Fig. [S] assuming 
that the FG of NGC 63 97 has been fully lost, as suggested by 
\D'Antona & Caloi hOOA) . 
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Figure 9. Results for the model NGC 2808-0. Data for giants 
(open circles) are taken from Carretta et al. (2009a, 2009b). The 
two turnoff stars examined by Bragagliaet al. (2010) are shown 
in the Mg— Al panel with star symbols. The open squares in the 
O-Na panel represent the location of stars formed from ejecta 
with Y>0.34 in absence of deep-mixing (e.g. for the blue MS 
turnoff stars) . The remaining symbols are as in Fig. |2] The top 
right panel shows the histogram of the helium distribution, with 
relative number fractions in evidence. 

other hand, require a larger star formation efficiencjQ. Mod- 
els with a short SG formation phase can be compatible with a 
second phase of SN II explosions due to the evolution of SG 
massive stars, thus relaxing any assumption of anomalous 
IMF. 

(iii) As an aside of the previous point, we stress that a 
shorter SG formation implies that the mass range of the 
AGB progenitors providing polluted ejecta for the SG forma- 
tion is narrower and, therefore, that these models require a 
more massive FG cluster to be able to produce the observed 
amount of SG stars. Models with a short SG formation phase 
typically require an initial FG mass 3-4 times larger than 
that needed in models with a more extended SG formation 
phase. 



5 RESULTS: NGC 2808 

We model NGC 2808 by taking it as the prototype of clusters 
with an extended Na-0 anticorrelations, along the lines de- 
scribed by D'Ercole et al. (2008) and in D2010. This means 
that we do not attempt to model the recent observational 
details in the analysis of the red HB stars, shown by Gratton 



Figure 10. As Fig.|9] but for model NGC 2808-1. 



The hydrodynamical simulations by iD'Ercole et al.l ||2008|) 
proved to be independent of the efficiency of star formation — 
see their Appendix A. This remains true when the timescale of 
star formation is long enough that all the gas is consumed. For 
models with shorter SG formation timescale the role of u becomes 
critical. 
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Figure 11. Results for the model NGC 2808-2 . The parameters 
are the same as to those of model NGC 2808-1 in Fig. 1101 but the 
super-AGB ejecta are assumed to accumulate without forming 
stars until t = 50 Myr, as can be seen from the bottom right 
panel (symbols as in Fig. [9]l. 



et al. (2011), and consider only the formation of two main 
populations: an extreme one, born from super-AGB undi- 
luted ejecta, and a milder one, born after dilution with the 
re-accreted pristine gas. 

5.1 Continuous star formation 

Figure [g] shows the model NGC 2808-0 which has a set of 
parameters very similar to that for the model of NGC 2808 
illustrated by D2010. After the phase of formation from 
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super-AGB ejecta, the peak value of re-accretion occurs at 
tacc=58 Myr (cf. the bottom right panel). The epoch from 
38.9 to ^ ~50 Myr corresponds to the formation of the 
Extreme population, with Y~0.35. This accounts for the 
presence of the blue main sequence. The following dilution 
phase leads to the formation of a less extreme SG (Y~0.3) 
that populates the intermediate main sequence of this clus- 
ter. 

The top left panel of figure [5] shows the Mg and Al 
predictions. Using the new yields for Mg and Al listed in 
Table 1 does not allow the fit of the three Mg-poor gi- 
ants f5fMg/Fel — 0.6) obser v ed by Carretta et al. 2009b. 
I Ventura. Carini &: D'Antonal (|201ll ') deal extensively with 
the Mg-Al problem in the models, and show that small ad- 
justments to the cross sections (concerning the Mg isotopes) 
may well reproduce t he [Mg/Fe]-[A1/Fe] va lue of the blue 
MS stars observed bv lBragaglia et al.l l|201(J ) and shown in 
Fig. [5] with the star symbol. In particular, an increase by 
a factor two of the cross section ^^Mg(p,7)^®Al shifts the 
6 Mq yield in the region occupied by the M 13 giants and 
by the MS star of NGC 2808. Such an increase is also con- 
sistent with recent new measurem ents of the reaction rate 
IlUNA Collaboration et al.ll201ll ). 

The difference in Mg between the MS stars and 
the giants can not be attributed to further dilution in 
the red giants, as magnesium is not touched in the 
red giant ev olution of low mas s stars, not even in 
the interior ([D'Antona fc Ventural [2007). Figure 5 in 
IVentura. Carini fc D'Antonal (|201ll ') shows the Mg-Al data 
in the literature for the clusters NGC 6752, M 15, NGC 
2808 and M 13, and the results of different modeling of the 
yields in the chemistry evolution of a 6 Mq star as adopted 
in this paper. Apart from two giants in M 15, a cluster whose 
metal abundance is much smaller than the others, no star 
has Mg abundances as low as those in the three UVES giants 
of NGC 2808, that are not compatible with the models. As 
these three stars constitute a problem for any other pollution 
model, and can not be explained in the deep extra mixing 
framework as well, we defer the study of this problem to a 
future investigation. 

Let us now consider the [Na/Fc]-[0/Fe] diagram (bot- 
tom left panel of Fig.[§]). While in the case of M 4-0 the new 
simulation is very similar to that presented in D2010 for the 
absence of an Extreme population, here the new super-AGB 
sodium yields are able to account for the high sodium values 
of the stars with upper limit only available for oxygen (full 
squares); this results ensues from the deep- mixing hypotesis, 
and the subsequent reduction of the oxygen values given by 
the recipe described in Sect. 12.21 The open squares represent 
the values of oxygen that should be found in turnoff stars for 
the stars with Y>0.34, as in these stars deep-mixing has not 
yet taken place. These stars, as well as blue MS stars, should 
be characterized by oxygen abundances larger than those of 
giant stars and fall along the Na-0 curve determined by the 
new yields and shown in Fig. [l] a prediction that can be 
falsified by spectroscopic observations. 

The top right panel of Fig. [5] shows the histogram of 
the distribution of helium abundances. The number ratio of 
the Intermediate and Extreme populations of the SG de- 
pends on the initial mass function and on the timing and 
extent of the dilution with pristine gas, and thus this is a 
powerful constraint for the successful models. The observed 



number ratios of the three populations differing in helium 
content in NGC 2808 is constrained by two different sets of 
observations: 



(i) the number counts of the HB stars (jBedin et al.ll200Gl ). 
These numbers were the basis of the interpretation of the HB 
in ter ms of multiple populations with different helium con- 
tent JD'Antona fc Caloi"2004': 'P'Ant ona et al.l l2005'l. pro- 
viding fractions of 50%, 35% and 15%, respectively, for the 
groups of standard helium (red clump stars), intermediate 
helium (blue luminous HB, named EBTl) and high helium 
(sum of the EBT2 and EBT3 extensions of the blue HB). 

(ii) t he number counts in the tripl e MS. [M ilonc ct ajj 
(|2011bl ') recently re-analyzed the data bv lPiottoet al.. (2007l 'l 
in order to derive the mass functions of the three popula- 
tions. According to this analysis, the fractions of standard 
(red MS), intermediate (middle MS) and high helium (blue 
MS) stars are respectively 62%, 24% and 14%. Therefore, 
while the fraction of FG stars changes from ~50 to ~60%, 
the proportion of Intermediate and Extreme population are 
very similar in the two analyses: the Extreme population is 
~35% of the whole SG. 

In the present models for NGC 2808 we set x = 0.4-0.5, 
that is, we assume that the fraction of FG stars, with stan- 
dard helium content Y=0.24, is 60-50%, and set the other 
input parameters in order to obtain a reasonable result for 
the number ratio of the Extreme and Intermediate SG. 

From the above discussion it is clear that, contrary to 
the case of M 4, there is less freedom in the choice of the 
parameters because of the presence of the Extreme popula- 
tion and the need to reproduce the right proportion among 
the Extreme, the Intermediate and the Primordial popula- 
tions (cf the three peaks in the He histogram in Fig. 19]). 
Some degree of flexibility is still present though. Figure [10] 
shows another model for NGC 2808, whose input is model 
NGC 2808-1 in Table 2. In the left panels we see that, as for 
the model NGC 2808-0, when dilution with pristine matter 
sets in, the model trajectories (dashed line) tend towards the 
FG abundances, but, during the last 10-15 Myr, star forma- 
tion occurs again in mostly undiluted AGB ejecta, and the 
abundances revert back to lower O and larger Na values. In 
fact, this last phase of star formation is necessary to fill the 
gap in O, Na values between the formation of the Extreme 
population (with very large Y, corresponding to the blue 
MS) and the Intermediate population with strong dilution. 
The bump in the helium distribution between the peaks at 
Y=0.3 and Y=0.35 (top right panel) forms during this same 
stage as well. And also the stars with very large [Al/Fe] in 
the simulation (upper left panel) form in this last epoch of 
star formation. This can be understood by looking at the 
abundances listed in Table 1, showing that the maximum 
abundance in Al and maximum depletion in Mg are reached 
with the evolution of the 5.5M0(at about 85 Myr). 



5.2 Star formation after gas accumulation 

We now focus on Fig.[TT]that shows the results of the model 
NGC 2808-2, having the same input parameters of model 
NGC 2808-1, but dealing with a different hypothesis for 
the star formation: here we assume that SG formation does 
not start immediately after the end of the SN II epoch, but 
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Figure 12. Case NGC 2808-3. Star formation is interrupted for 
30 Myr, but tlie second stage of star formation lasts for 20 Myr 
(symbols as in Fig.[9]l. 



that gas accumulates from 40 to 50 Myr, and then a star 
formation burst occurs. 

The results are very similar, but the Extreme popula- 
tion does not reach the high sodium abundances obtained 
in the model NGC 2808-1; this is a consequence of the fact 
that, before the SG formation starts, the ejecta of the stars 
between 8 Mq and ~7.3 MquiIx, so that the highest abun- 
dances are lowered. 



5.3 Models with interruption of star formation 
due to SG SN II explosions 

As discussed in Sect.U an interesting implication of clusters 
like M 4, with a short [0/Fe]-[Na/Fe] anticorrelation, is that 
thay can form the bulk of the SG population in 10 Myr, or 
even less; such a timescale is compatible with the survival of 
a cooling flow in the cluster core, previous to the explosion 
of the Type II SN o/ the second generaUoqj, and thus there 
would be no reason to limit the IMF of the SG to the stars 
that do not explode as SN II. 

Understanding whether the formation of two distinct 
SG populations (such as those found in NGC 2808) is com- 
patible with star formation bursts lasting ~10 Myr is not 
straightforward. We have tried to model in detail the three 
populations of NGC 2808 by making the assumption that, 
~10 Myr after the Extreme population began to be formed 
in the cluster, a new phase of SN II begins. This will hamper 
star formation until the end of the SN II phase. 

In our models, this phase lasts for ^ 38 Myr (Table 1). 



^ Although the most massive stars may take only 5 Myr to ac- 
complish their complete evolution, it is possible that only stars 
below a mass evolving in ^10 Myr produce a SN II explosion, 
the more massive stars may end their life in the less energetic 
formation of a black hole llSmarttil2009l ) 
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Figure 13. Case NGC 2808-4 (top) and 5 (bottom). Using the 
input parameters of Fig. 1121 but limiting the second star forma- 
tion phase to 10 Myr, the Intermediate population is very scarce 
(top panel). Increasing the infall of pristine matter, the O-Na 
data come very close to the FG values (bottom panel) (symbols 
as in Fig.|9ll. 



We can have a second burst of star formation from the AGB 
ejecta starting at ^^90 Myr. We have discussed in Sect.[2]that 
the time hiatus between the two bursts may be reduced to 
~30 Myr for models employing a smaller core-overshooting. 
Thus we employ in the following ~30 Myr as total duration 
of the SG SN II epocljl. 



Notice that the timing of each mass for our simulations, given 
in Table 1, is still adopted in the same way in these models. 
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The main difficulty encountered in modeis inciuding 
SN II expiosions from the Extreme SG is to reproduce at 
the same time the extension and distribution of stars aiong 
the O-Na anticorreiation and the heiium distribution func- 
tion. 

• Case NGC 2808-3 (Fig. [H)): In this model we stop 
the SG formation from pure ejecta at 50 Myr, allow for a 

^ 30 Myr period without star formation, and then resume 
star formation from the AGB ejecta mixed with pristine gas, 
with a choice of tacc=87 Myr. The clumps in the distribution 
of stars in the O-Na plane are probably too separated to be 
consistent with the observations. Moreover, in order to ob- 
tain a reasonable number ratio between the Extreme stars 
(Y>0.34) and the Intermediate population (Y~0.30), we are 
forced to extend the star formation of this latter population 
to ~20 Myr, much longer than allowed by the beginning of 
the new (third) SN II epoch. 

• Case NGC 2808-4 (Fig. [l3 top panel): this model is 
identical to NGC 2808-3, but we limit the formation of the 
Intermediate population to 90 Myr: as discussed above, the 
Intermediate population is a bare 11%, too small to repro- 
duce the number ratio of the Extreme- and Intermediate- 
heUum MS. 

• (Case NGC 2808-5, Fig. 1131 bottom panel): we increase 
the relative number of Intermediate population stars by in- 
creasing the dilution. The second peak of helium content 
occurs at Y^O.28, and the number ratio of Extreme and 
Intermediate population becomes reasonable, but the O- 
Na values become much too close to the F G values, due to 
the st rong dilution (a problem discussed in iD'Ercole et al.l 
l|201ll ) in a different context). 

We find then that these simulations do not reproduce 
the observations so well as models NGC 2808-1 and 2. 
Nevertheless, the variety of observational results for dif- 
ferent clusters may have their explanation in the variety 
of relevant parameters. For instance, model NGC 2808-4 
(Fig. [T3] top panel) could be relevant to explain the hori- 
zontal branch stellar dist ribution of the cluster NGC 2419 
l|di Criscienzo et al.|[201ll '). 



6 CONCLUSIONS 

In this paper we have explored the possible role of super- 
AGB in the formation of the observed abundance patterns 
in multiple-population globular clus ters. Our models are 
based on the set of yields calculated bv lVentura fc D'Antonal 
((20091, I2OIII ) . The results of our investigation show that the 
new super- AGB yields broaden the range of possible forma- 
tion histories leading to the observed abundance patterns. 
Specifically: 

• we have built models for M 4, a cluster representative 
of numerous other GCs characterized by a shorter Na-O 
anticorreiation. For this cluster, we show that the observed 
abundance pattern can be reproduced also by models with 
a short SG formation phase ('^ 10-20 Myr) involving only 
super-AGB ejecta and pristine gas. We confirm that models 
presented in our previous work (D'Ercole et al. 2010) which 
are characterized by a more extended SG formation epoch 
involving super-AGB and AGB ejecta (along with pristine 
gas) can also reproduce the observed abundance patterns. 



While models with a shorter SG formation epoch can in- 
clude the effect of SG SN II, the narrower range of AGB 
progenitors providing gas for SG formation implies that they 
require a larger FG initial mass (by a factor 3-4 compared 
to models with a more extended SG formation phase). 

• In order to account for the O-poor Na-rich stars of NGC 
2808 we assumed the presence of deep-mixing in SG giants 
forming in a gas with helium abundance Y> 0.34, which 
significantly reduces the atmospheric oxygen content, while 
preserving the sodium abundance. The use of the super- 
AGB yields does not change the main requirements on SG 
formation history of models for the more complex clusters 
like NGC 2808. In these clusters, an Extreme population 
with very large helium content form from the pure ejecta of 
super-AGB stars, followed by formation of an Intermediate 
population by dilution of massive AGB ejecta with pristine 
gas. We also attempted to model the three populations of 
NGC 2808 allowing for a long hiatus between the formation 
of the Extreme and Intermediate population (possibly due 
to SN II explosions of the Extreme populations), but our 
attempts to model each star formation epoch in a period as 
short as ~10 Myr were not satisfactory. Nevertheless, further 
study is probably needed to settle this issue. 

We finally want to point out that one-zone models such 
as those presented here (and in D2010), while very useful to 
understand rather easily and very quickly the role of differ- 
ent ingredients in shaping the GC chemical evolution (AGB 
ejecta, amount of pristine gas, duration of the system evo- 
lution, etc.), cannot of course capture the complexity and 
variety of effects connected with the gas hydrodynamics. As 
a simple example, we mentio n that, during the co oling flow 
evolution of the GC gas (c.f. ID'Ercole et aklboOSl ). most of 
it cumulates in the centre of the cluster which therefore at 
each time hosts not only the AGB ejecta delivered "in situ" , 
but also that delivered at earlier times by FG stars located 
at the GC outskirts. It is thus expected that SG stars with 
different chemical properties can form in different places at 
the same time, and not only at different times as in one-zone 
models. For this reason w e plan to work out hy drodynamic 
models similar to those of lD'Ercole et al.l ([2003), but taking 
into account the detailed chemistry as in the present models. 
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